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Abstract 

The azimuthal distribution of jets produced in the Breit frame in high-Q^ deep 
inelastic e'^p scattering has been studied with the ZEUS detector at HERA using 
an integrated luminosity of 38.6 pb~^. The measured azimuthal distribution 
shows a structure that is well described by next-to-leading-order QCD predictions 
over the range considered, > 125 GeV^. 
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1 Introduction 



Jet production in neutral current (NC) deep inelastic scattering (DIS) at high (where 
is the negative of the virtuality of the exchanged boson) provides a testing ground 
for the theory of the strong interaction between quarks and gluons, namely quantum 
chromodynamics (QCD). An observable of interest is </)j^^, the azimuthal angle in the 
Breit frame [1] between the lepton scattering plane, defined by the incoming and outgoing 
lepton momenta, and the jets produced with high transverse energy {E^-^^^) in that frame. 

In the Standard Model, azimuthal asymmetries arising from perturbative QCD effects [2, 3, 4] 
are expected in the distribution. At leading order (LO), the azimuthal dependence 
for unpolarised NC DIS at <^ M| has the form 

^ = A + 5 cos (0j^J + C cos (20jj . (1) 

The current-current form of the electromagnetic interactions makes the cross section linear 
in cos(0j^J, cos(2</)j^J, sin(</)j^J and sin(2{/)j^J. However, the coefficients of the terms in 
sin((/)j^^.) and sin{2(j)^^^) vanish due to time-reversal invariance and the absence of final- 
state interactions at the quark-gluon level at LO. The coefficients A, B and C result 
from the convolution of the matrix elements for the partonic processes with the parton 
distribution functions (PDFs) of the proton [3,4]. The cos(2(/)j^^) term is expected from 
the interference of amplitudes arising from the -|-1 and —1 helicity components of the 
transversely polarised part of the exchanged photon, whereas the interference between 
the transverse and longitudinal components gives rise to the cos(</)j^^-) term. In addition, 
a non-perturbative contribution to the asymmetry arises from the intrinsic transverse 
momentum of partons in the proton. Since such intrinsic transverse momenta are small [5], 
this contribution is expected to be negligible for jet production at high -E-rjct [6]- 

Previous studies of single hadron production in NC DIS observed a cos0 term that was 
attributed to non-perturbative effects [7]. However, more recently, a ZEUS measure- 
ment of the azimuthal dependence of charged hadrons with high transverse momentum in 
the centre-of-mass system gave evidence for perturbative contributions to the azimuthal 
asymmetry [8]. This paper presents the first study of the azimuthal distribution of jets 
with high transverse energy in the Breit frame and the comparison with LO and next-to- 
leading-order (NLO) QCD predictions. 



2 Experimental details 

These results are based on data collected in 1996-1997 with the ZEUS detector at HERA, 
corresponding to an integrated luminosity of 38.6 ± 0.6 pb~^. The HERA rings were op- 
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erated with protons of energy Ep = 820 GeV and positrons of energy = 27.5 GeV. The 
ZEUS detector is described elsewhere [9,10]. The main components used in the present 
analysis are the central tracking detector [11], positioned in a 1.43 T solenoidal magnetic 
field, and the uranium-scintillator sampling calorimeter (CAL) [12]. The tracking detec- 
tor was used to establish an interaction vertex. The CAL covers 99.7% of the total solid 
angle. It is divided into three parts with a corresponding division in the polar angleQ, 6, 
as viewed from the nominal interaction point: forward (FCAL, 2.6° < 9 < 36.7°), barrel 
(BCAL, 36.7° <9 < 129.1°), and rear (RCAL, 129.1° <9 < 176.2°). The smallest subdi- 
vision of the CAL is called a cell. Under test-beam conditions, the CAL relative energy 
resolution is 18%/ ^^/E{GeV) for electrons and 35%/ ^y E{GeV) for hadrons. A three- level 
trigger was used to select the events online [10]. 

As the analysis follows very closely that of the inclusive jet cross sections in the Breit 
frame [13], details about the event selection, jet finding, systematic uncertainties and 
theoretical predictions are not repeated here. 

The scattered-positron candidate was identified from the pattern of energy deposits in 
the CAL [14]. The kinematic region of the analysis was selected by the requirements 
> 125 GeV^ and — 0.7 < cos 7 < 0.5, where 7 is the angle of the scattered quark in 
the quark-parton model. Cuts on this angle restrict the phase-space selection in Bjorken 
X and the inelasticity y due to the relation 

(1 - y)xEp - yEe 

cos 7 = 7 r — — . 

(1 - y)xEp + yE^ 

The longitudinally invariant cluster algorithm [15] was used in the inclusive mode [16] 
to reconstruct the jets in the hadronic final state both in data and in events simulated by 
Monte Carlo (MC) techniques. In data, the algorithm was applied to the energy deposits 
measured in the CAL cells after excluding those associated with the scattered-positron 
candidate. The jet search was performed in the pseudorapidity (?7'^)-azimuth {4>^) plane 
of the Breit frame, where (1)^ = corresponds to the direction of the scattered positron. 
The transverse energy of the jets in the Breit frame, E^^^^, was required to be larger than 
8 GeV and the pseudorapidity range was restricted to —2 < r]^^ < 1.8. The data sample 
contained 8523 events, 5073 of which were one-jet, 3262 two-jet, 182 three-jet and 6 four- 
jet events. The range covered by the data sample extended up to ~ 16 000 GeV^; 
measurements of the azimuthal distribution are presented up to a mean value of 
~ 2300 GeV^. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is 
defined as 77 = — ln(tan |), where the polar angle 9 is taken with respect to the proton beam direction. 
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3 Monte Carlo studies and systematic uncertainties 



Samples of events were generated to determine the response of the detector to jets of 
hadrons and to calculate the correction factors necessary to obtain the hadron-level jet 
cross sections. The generated events were passed through the GEANT 3.13-based [17] 
ZEUS detector- and trigger-simulation programs [10] and were reconstructed and anal- 
ysed by the same program chain as the data. The NC DIS events were generated using 
the LEPTO 6.5 program [18] interfaced to HERACLES 4.6.1 [19] via DJANGOH 1.1 [20]. 
The HERACLES program includes photon and Z exchanges and first-order electroweak 
radiative corrections. The QCD cascade was modelled with the ARIADNE 4.08 pro- 
gram [21]. The CTEQ4D [22] parameterisations of the proton PDFs were used. As an 
alternative, samples of events were generated using the model of LEPTO based on first- 
order QCD matrix elements plus parton showers (MEPS). In both cases, fragmentation 
into hadrons was performed using the JETSET 7.4 program [23]. In these programs, the 
azimuthal distribution was generated according to the LO QCD calculation. 

The comparison of the reconstructed jet variables for the hadronic and the calorimetric jets 
in simulated events showed that no correction was necessary for and that the average 
resolution was 0.09 radians. The sample of events generated with either ARIADNE or 
LEPTO-MEPS, after applying the same offline selection as for the data, gave a good 
description of the measured distributions for both the event and jet variables [13,24]. 
However, a MC sample of events generated with a uniform azimuthal distribution did not 
describe the observed distribution at detector level. These comparisons establish the 
presence of an azimuthal modulation in the data. 

The cross sections presented here were corrected to the hadron level by applying bin-by-bin 
corrections to the measured distributions. The correction factors had some dependence 
on due to the cuts applied to remove the effects of QED radiation that lead to a 
radiated photon from the positron being misidentified as a hadronic jet. The observed 
dependence of the correction factor was not sensitive to the assumed azimuthal distri- 
bution in the generator; this was confirmed by the observation that the correction factors 
based on the MC sample generated with a uniform azimuthal distribution had the same 
dependence on The MC programs were also used to evaluate the corrections for QED 
radiative effects, which were negligible for the normalised cross sections. 

A detailed study of the systematic uncertainties was carried out. Those that had an effect 
on the shape of the azimuthal distribution were: 

• the uncertainty in the absolute energy scale of the jets; 

• the uncertainty in the MC modelling of the hadronic final state, which was estimated 
from the differences between ARIADNE and LEPTO-MEPS in correcting the data for 
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detector effects; 



• the uncertainty in the positron identification, which was estimated by repeating the 
analysis using an alternative technique [25] to select the scattered-positron candidate. 

The relative changes in the normalised differential cross section induced by the variations 
mentioned above were typically smaller than the statistical uncertainties, which ranged 
from ~ 2% at _ 125 GeV^ up to 6% at ^ iqoo GeV^. 



4 Perturbative QCD calculations 



The LO and NLO QCD predictions were obtained using the program DISENT [26]. 
The number of flavours was set to five and the renormalisation {fin) and factorisation 
(/ii?) scales were chosen to be fi^ = E^^-^t I^f = respectively. The strong cou- 
pling constant, a^, "was calculated at two loops with A-^ = 220 MeV, corresponding to 
as{Mz) = 0.1175. The calculations were performed using the MRST99 [27] parameter- 
isations of the proton PDFs. The results obtained with DISENT were cross-checked by 
using the program D1SASTER++ [28]. The differences were always smaller than 1%. 

The perturbative QCD contribution to the terms B and C in Eq. (|Tp is large. At LO 
in as, two processes contribute to jet production in the Breit frame: QCD-Compton 
scattering (QCDC, 7*g qg) and photon-gluon fusion (PGF, 7*(7 qq). For the 
former, the scattered gluon (quark) preferentially appears at ^j^^ = (vr), whilst for the 
PGF process, the (j)^^^ dependence is dominated by the cos(20j^^.) term and is very similar 
for quarks and antiquarks. Thus, the different contributions to the cos{(j)f^^ term from 
quarks and gluons tend to cancel in the cos(0j^^) asymmetry and the predicted azimuthal 
distribution is very close io A + C cos(2</)j^^). The NLO QCD correction mainly modifies 
the normalisation and slightly affects the shape of this prediction. In order to test the 
QCD prediction for the azimuthal distribution, it is desirable that no cut be applied to 
the jets in the laboratory frame; otherwise, the azimuthal distribution can be strongly 
distorted by kinematic effects [4]. For this reason, no such cut was used in the definition 
of the cross sections. 

Since the measurements refer to jets of hadrons, whereas the perturbative QCD calcula- 
tions refer to partons, the hadronisation effects were investigated by using the models of 
ARIADNE, LEPTO-MEPS and HERWIG [29]. These effects were negligible [24]. 
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5 Results 



The cross sections presented here include every jet of hadrons in an event with -E-rjct > 
8 GeV and —2 < 7]^^^ < 1.8. A detailed comparison of the differential cross sections as 
functions of Q^, E^-^ot ^"^^ V^t inclusive jet production in the same kinematic region as 
used here was presented in a previous publication [13]. At low and low E^jct^ NLO 
QCD calculations fall below the data by ~ 10%. Nonetheless, the differences between the 
measurements and calculations are of the same size as the theoretical uncertainties [13]. 
The comparison of the shape of interest in this publication is facilitated by normalising 
the predicted cross section and the data to unity. 

The normalised differential cross-section (l/cr) da/d(f)f^^ for inclusive jet production as a 
function of (f)f^^ is shown in Fig. |l| and in Table |l[ This distribution has clear enhancements 
at = and ^j^^ = vr. The NLO QCD calculations with either /i/j = E^^^^ or Q 
reproduce the asymmetry. This comparison constitutes a precise test of the perturbative 
QCD prediction for the azimuthal distribution since the theoretical uncertainties are small. 
The dominant theoretical uncertainty arose from terms beyond NLO and was estimated 
by varying /ir between E^-^^^/2 and 2E^-^^^; the effect on the amplitude of the modulation 
of the distribution was ~ ±1%. This observation complements the ZEUS measurement 
of the azimuthal dependence of charged hadrons with high transverse momentum in NC 
DIS [8]. 

The measurements folded about vr, in different regions of are presented in Fig. ^ 

and in Table ^ The LO and NLO QCD predictions are compared to the data. The 
NLO QCD prediction describes the data well, whereas the LO QCD calculations predict 
a larger asymmetry, particularly in the lower intervals. In both cases, the asymmetry 
is predicted to decrease as increases, as a result of the progressive decline of the 
contribution from the PGF process. 

In order to perform a more quantitative study of the asymmetry and its dependence on 
Q^, a fit was performed to the values of (l/o") da / d\(l)^^\ both in the data and in the QCD 
predictions. The functional form 

l + /icos(0j^j + /2cos(20j^j 

was used. The parameter fi (/2) represents the contribution of the cos^j^^ (cos20j^J 
term to the total asymmetry. The fitted values of /i and /2 as functions of and for 
the entire sample with > 125 GeV^ are shown in Fig. |] and listed in Table together 
with the LO and NLO QCD predictions and their uncertainties. The comparison of the 
LO QCD calculations for the QCDC and PGF process shows that the asymmetry is 
predicted to arise mostly from the gluon-induced interactions. The LO QCD predictions 



1 da _1 
a d\(l)fj 71 
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do not reproduce the measurements. However, the uncertainty at LO is rather large. The 
difference between the LO and NLO calculations has been assigned as the theoretical 
uncertainty of the LO predictions and is ~ ±0.04 (±0.01) for /2 (/i). At NLO, the 
dominant theoretical uncertainty on /2 (/i) was that due to terms beyond NLO and was 
estimated by fitting the predictions obtained with fi^ = E!^^^^^/2 and 2El^-^^{, it amounted 
to ~ ±0.01 (±0.005). The NLO predictions for /i and /2 based on calculations using 
f^R = Q differed from those using fi^ = -E-rjct t>y as much as the estimated theoretical 
uncertainty. The NLO QCD predictions are in good agreement with the measured values 
of f2- For /i, the observed asymmetry tends to be slightly larger and more negative 
than that predicted by perturbative QCD. The measurements are consistent with the 
dependence of /i and /2 predicted by NLO QCD. 



6 Summary 

A study of the azimuthal asymmetry for inclusive jet production in neutral current deep 
inelastic e~^p scattering in the Breit frame at a centre-of-mass energy of 300 GeV has been 
presented. Jets of hadrons were identified with the longitudinally invariant kx cluster 
algorithm in the Breit frame. The normalised cross sections as a function of the azimuthal 
angle of the jets in the Breit frame are given in the kinematic region > 125 GeV^ and 
—0.7 < cos 7 < 0.5. The cross sections include every jet of hadrons in the event with 
El^j^t > 8 GeV and -2 < r^f^, < 1.8. 

The measured azimuthal distribution peaks in the directions along, and opposite to, that 
of the scattered positron in the Breit frame. The NLO QCD calculations give a good 
description of the observed azimuthal variation. The dependence of the azimuthal asym- 
metry on is also compatible with NLO QCD. 

These measurements constitute a precise test of the perturbative QCD prediction for the 
azimuthal distribution since the theoretical uncertainties are small. 
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interval (rad) 




rirr /rlA^^ 

"^/"<Pjct 


A 


A 

'—^syst 


0- 


7r/6 





1655 


±0.0054 


+0.0042 
-0.0015 


7r/6 — 


7r/3 





1630 


±0.0051 


+0.0011 
-0.0014 


7r/3- 


7r/2 





1398 


±0.0047 


+0.0020 
-0.0008 


7r/2 - 


27r/3 





1557 


±0.0050 


+0.0000 
-0.0022 


27r/3- 


57r/6 





1601 


±0.0051 


+0.0062 
-0.0013 


57r/6 — 


TT 





1771 


±0.0057 


+0.0050 
-0.0075 


TT — 


77r/6 





1779 


±0.0056 


+0.0016 
-0.0052 


77r/6- 


47r/3 





1577 


±0.0050 


+0.0051 
-0.0015 


47r/3- 


37r/2 





1458 


±0.0046 


+0.0035 
-0.0008 


37r/2- 


57r/3 





1468 


±0.0047 


+0.0032 
-0.0028 


57r/3- 


ll7r/6 





1531 


±0.0047 


+0.0017 
-0.0014 


ll7r/6- 


27r 





1674 


±0.0054 


+0.0018 
-0.0027 



Table 1: Normalised differential cross-section (1/cr) da/dcp^^ for inclusive jet 
production with -ETjet > 8 GeV and —2 < 77-^^ < 1.8. The statistical and systematic 
uncertainties are shown separately. 
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interval 



(rad) 



{I /a) da/d\d^f,,\ 



A 



stat 



A 



syst 



{I /a) da/d\cj,f,,\ 



A 



stat 



A 



syst 



125 < g2 < 250 GeV^ 



250 < Q2 < 500 GeV^ 



7r/6 
7r/3 





7r/6 
7r/3 -7r/2 
7r/2 - 27r/3 

27r/3 - 57r/6 

57r/6 



TT 



0.3319 
0.3171 
0.2932 
0.2907 
0.3232 
0.3538 



±0.0103 
±0.0096 
±0.0095 
±0.0093 
±0.0101 
±0.0109 



+0.0068 
-0.0069 

+0.0054 
-0.0028 

+0.0066 
-0.0085 

+0.0038 
-0.0018 

+0.0120 
-0.0021 

+0.0049 
-0.0094 



3461 
3116 
2754 
3259 
3011 
3497 



±0.0138 
±0.0122 
±0.0116 
±0.0126 
±0.0122 
±0.0141 



+0.0071 
-0.0073 

+0.0054 
-0.0072 

+0.0060 
-0.0047 

+0.0032 
-0.0051 

+0.0129 
-0.0033 

+0.0059 
-0.0126 



500 < g2 < 1000 GeV^ 



g2 > 1000 GeV^ 



- 7r/6 
7r/6 — 7r/3 
7r/3 -7r/2 
7r/2 - 27r/3 

27r/3 -57r/6 

57r/6 — vr 



0.3268 
0.3136 
0.2713 
0.2871 
0.3418 
0.3693 



±0.0192 
±0.0178 
±0.0164 
±0.0167 
±0.0187 
±0.0206 



+0.0100 
-0.0085 

+0.0063 
-0.0055 

+0.0053 
-0.0052 

+0.0079 
-0.0062 

+0.0075 
-0.0036 

+0.0062 
-0.0120 



3129 
3210 
3126 
2989 
3074 
3571 



±0.0229 
±0.0220 
±0.0211 
±0.0202 
±0.0215 
±0.0247 



+0.0064 
-0.0047 

+0.0068 
-0.0182 

+0.0177 
-0.0039 

+0.0027 
-0.0009 

+0.0178 
-0.0048 

+0.0105 
-0.0299 



g2 > 125 GeV^ 



- 7r/6 
7r/6 — 7r/3 
7r/3 -7r/2 
Tx/2 - 27r/3 

27r/3 - 57r/6 

57r/6 



TT 



0.3334 
0.3153 
0.2867 
0.3016 
0.3176 
0.3552 



±0.0072 
±0.0066 
±0.0064 
±0.0065 
±0.0068 
±0.0076 



+0.0053 
-0.0043 

+0.0026 
-0.0027 

+0.0035 
-0.0019 

+0.0017 
-0.0014 

+0.0116 
-0.0025 

+0.0044 
-0.0117 



Table 2: Folded normalised differential 
regions of for inclusive jet production 
For details, see the caption of Table 



cross-section {1 / a) da / d\(j)^^^\ in different 



with E^.^, > 



8 GeV and 



-2 < r/^, < 1.^ 



10 





Q"^ region (GeV^) 




^stat 


^syst 


(PGF,QCDC) 


NLO QCD 


fl 


All Q2 (g2 > 125) 


-0.0273 


±0.0144 


+0.0121 

—0 ooqq 


0.0115 ±0.0118 
(0.0236,-0.0013) 


-0.0003 




125 < Q2 < 250 


-0.0248 


±0.0208 


+0.0113 
— u.uuyo 


0.0171 ±0.0100 

(0.0303,-0.0005) 


0.0071 

— U.UUOiJ 




250 < Q2 < 500 


-0.0103 


±0.0268 


+0.0144 

— VJ.UIDO 


0.0106 ±0.0136 
(0.0210,-0.0015) 


-0.0030 




500 < Q2 < 1000 


-0.0690 


±0.0388 


+0.0166 
— U.UioU 


0.0060 ±0.0161 

(0.0152,-0.0029) 


-0.0101 

— U.UUD ( 




Q2 > 1000 


-0.0238 


±0.0465 


+0.0196 
—0.0168 


0.0022 ±0.0122 

(0.0089,-0.0009) 


-0.0100 




All Q2 (q2 > 125) 


0.0947 


±0.0143 


+0.0068 
-0.0133 


0.1340 ±0.0356 
(0.1999,0.0452) 


0.0984 1°:°?^^ 




125 < Q2 < 250 


0.0969 


±0.0207 


+0.0095 
-0.0151 


0.1418 ±0.0388 
(0.1880,0.0410) 


1 030 +00074 

U.iUOU _o.0127 




250 < Q2 < 500 


0.0906 


±0.0270 


+0.0112 
-0.0164 


0.1496 ±0.0424 
(0.2262,0.0632) 


n 1 nvo +0.0088 

U.IU/Z _o.oi58 




500 < Q2 < 1000 


0.1348 


±0.0374 


+0.0044 
-0.0082 


0.1306 ±0.0356 
(0.1982,0.0358) 


09^0 +° °°^^ 

U.UJOU _o.oi46 




Q2 > 1000 


0.0526 


±0.0462 


+0.0086 
-0.0387 


0.0754 ±0.0160 
(0.1678,0.0359) 


n n^Q4 +0-0041 

U.UOy^ -0.0076 



Table 3: Measured values of the parameters fi and f2 in the different 
regions. The LO and NLO QCD predicted values calculated using DISENT and 
the MRST99 parameterisation of the proton PDFs are shown for comparison. The 
quoted uncertainties in the theoretical predictions are described in the text. 
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CO .31. 



J 0.2 - • ZEUS 96-97 




0.1 h 

NLO QCD: 

DISENT MRST99 (^lR=E?Je,) 

DISENT MRST99 (^lR=Q) 

Q I I I I I I I I I I I I I I I I I I I I 

7C/2 7U 37C/2 2% 

*fe, (rad) 



0.05 - 



Figure 1: The normalised differential cross-section {\ / a) da / d())f^^ for inclusive jet 
production with E^-^ct > 8 GeF and —2 < rj^^^ < 1.8 (points). The inner error bars 
represent the statistical uncertainty. The outer error bars show the statistical and 
systematic uncertainties added in quadrature. The NLO QCD calculations using 
DISENT and the MRST99 parameterisations of the proton PDFs are shown for 
two choices of the renormalisation scale. 
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0.2 
0.4 



0.3- 



0.2- 



0.4- 



0.3- 



ZEUS 

I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I 

^ 0.4 1- 125<Q^<250GeV^ _|_ 250 < Q'' < 500 GeV 



I I I I I I I I I I I I I I l " l I I I I I I I I I I I I I I ' 



500 <Q^< 1000 GeV^ 




I I I I I I I I I I I I 



AIIQ' 




0.2 I-. ■ ■ I ■ ■ ■ I ■ ■ ■ I 



■ ■ 



qSiOOO GeV^ 




z - 1 I I I I I I I I I I I 



-l—L. 



|(^fet I (rad) 

• ZEUS 96-97 

NLO QCD DISENT (^R=Er%,) 

- - - LO QCD DISENT (^R=Er%,) 



7C/2 



7C 



W%x I (""ad) 



Figure 2: The folded normalised differential cross-section (1/cr) da/d\(f)^^^\ for 
inclusive jet production with -E^jet > 8 GeV and —2 < rj^^ < 1.8 in different 
Q"^ regions (points). The inner error bars represent the statistical uncertainty. 
The outer error bars show the statistical and systematic uncertainties added in 
quadrature. The LO and NLO QCD calculations using DISENT and the MRST99 
parameterisations of the proton PDFs are also shown. 
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ZEUS 



- a) 



0.1 - 



-0.1 - 



f 



2 0.2 - b) 



0.1 - 



All 



Asymmetry from the cos (|)?, term 
• ZEUS 96-97 

- - - LO QCD prediction (H[,=E?jet) 

M NLO QCD prediction (E?,jet/2 < Hr < 2E?jet) 



Asymmetry from the cos 2(|)?, term 
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J I L 
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10 



(GeV^) 



Figure 3: The fitted values of a) fi and b) f2 from the folded normalised differential 
cross-section (l/a) da/d\(j)f^^\ 



for inclusive jet production with -E'Tjet > ^ GeV and 



—2 < r^j^^ < 1.8 as a function ofQ^ (points). The fitted values for the entire sample, 
> 125 GeV^, are shown on the left-hand side. The inner error bars represent 
the statistical uncertainty. The outer error bars show the statistical and systematic 
uncertainties added in quadrature. The results of the fits to the LO and NLO QCD 
predictions using DISENT and the MRST99 parameterisations of the proton PDFs 
are shown. The shaded bands represent the uncertainty of the calculations due to 
the dependence on the renormalisation scale. 
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